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ABSTRACT 



Context. Studying the appearance and properties of bipolar winds is critical to understand the stellar evolution from the AGB to the 
planetary nebula (PN) phase. Many uncertainties exist regarding the presence and role of binary stellar systems, mainly due to the 
deficit of conclusive observational evidences. 

Aims. We investigate the extended equatorial distribution around the early bipolar planetary nebula M 2-9 ("Minkowski's Butterfly 
Nebula") to gather new information on the mechanism of the axial ejections. 

Methods. Interferometric millimeter observations of molecular emission provide the most comprehensive view of the equatorial mass 
distribution and kinematics in early PNe. Here we present subarcsecond angular-resolution observations of the 12 CO J = 2-1 line and 
continuum emission with the Plateau de Bure interferometer. 

Results. The data reveal two ring-shaped and eccentric structures at the equatorial basis of the two coaxial optical lobes. The two 
rings were formed during short mass-loss episodes (~ 40 yr), separated by ~ 500 yr. Their positional and dynamical imprints provide 
evidence of the presence of a binary stellar system at the center, which yields critical information on its orbital characteristics, 
including a mass estimate for the secondary of < 0.2 M Q . The presence of a stellar system with a modest-mass companion at the center 
of such an elongated bipolar PN strongly supports the binary-based models, because these are more easily able to explain the frequent 
axisymmetric ejections in PNe. 

Key words. (Stars:) circumstellar matter - Stars: AGB and post-AGB - Radio lines: stars - Stars: mass-loss 



1. Introduction 

The appearance and shaping of bipolar winds in a stage when 
stars transit from the late asymptotic giant branch (AGB) to 
the planetary nebula (PN) phase is one of the most intrigu- 
ing open questio ns in stellar evolution. Studying bipolar winds 
(Balick & Frank 2002) is critical for unraveling this late stel- 
lar transition and ultimately for understanding the formation 
of PNe, which o ften show strongly axisymmetric shapes (e.g. 
ISahai et aL I l2007h . One mechanism that is supposed to trigger 
the generation of bipolar winds is the presence of a binary stellar 
system at the nebulae' s center (ISokerll200lb iFrank & Blackman 
20041 iDe Marcol l2"009). The direct detection of stellar compan- 
ions remains an observational challenge, however. Indeed, cen- 
tral stars of bipolar post-AG B nebulae are rarely found to be 
multiple (iHrivnaket al.ll2011l) . 

M 2-9 (also known as "Minkowski's Butterfly Nebula" or 
"the Twin Jet Nebula") is a prototypical young planetary neb- 
ula that displays a very elongated bipolar structure (120"xl2" 



in size, Schwa rz et al. I fl997l) . It became very popular from the 
spectacular images obtained by the Hubble Space Telescope 
(Bali ck et al.lll997t iBalickll 1999b in the central 60"-size region 
of the nebula, which show two coaxial th in outflows. In the equa- 
tor, the waist of the nebula is very slim dBalick et al. 1997), and 
its optical image is diffuse due to a dense torus of gas and dust. 

A sequence of obs ervations of ionized gas emission ob- 
tained every 2-5 years dDovle et al.ll2000t ICorradi et al.ll201 ll) 
shows a pattern of changes with mirror symmetry, with re- 
spect to the equatorial plane, that rotates around the symme- 
try axis of the nebula. A rotating ionizing/exciting beam from 
the star is postulated to explain the observed increase in bright- 
ness in a line of knots along the lobe edges. These observa- 
tions indicate that there is a binary stellar system at t he center 
dSchmeia & Kimeswengerf2001l;ISmith & Gehrzl2005l). orbiting 
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with a period ~ 90 years dDovle et al.ll2000tlCorradi et alj|2rjTu 
iLivio & SokerfeOOlUSmith et al.ll2005l) . " 

Several theoretical studies have addressed the possible rela- 
tionship between nebular bipolarity and binarity of its nucleus, 
which is one of the major long-standing problems toward a bet- 
ter understanding of stellar evolution. Particularly, for M 2-9 
it was proposed dLivio & Sokerll200T|) that a white dwarf star 
accretes mass from the wind of an AGB or post-AGB compan- 
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Fig. 1. Channel maps of the 12 CO 7 = 2-1 line em ission toward M 2 -9 (in contours) superimposed on the optical image obtained 
with the Hubble Space Telescope (in color scale; iBalick et alJI 1997h . The center is given by the position of compact continuum 
emission, here subtracted, at the J2000 coordinates RA 17:05:37.958, Dec -10:08:32.48. The LSR velocities are specified in the top- 
left corner of each panel. Contours are shown from 4 cr with a spacing of 5 cr (where the root-mean-square noise is 7.9 mJy beam" 1 ) 
in solid line, and at -4 cr in dashed lines. The synthesized beam is 0786 x 0740 at PA = 15° and is drawn in the bottom-right corner 
of the last panel. 



ion, and that a subsequent increase in magneto-centrifugal forces 
launches the fa st bipolar jets responsible for the rotating knots . 
Since the core dLvkou et alJl201 it iTorres-Peimbert et alJ feOlO) 
is embedded within a high-density dusty environment, there is 
no direct evidence of this scenario yet. On the other hand, the 
presence of disks in rotation in the innermost equatorial re- 
gions of se veral young PNe has been found to be associated 
to binarity dBuiarrabal et alj|2005b Ivan Winckelll2003l) . Several 
recent works have inves tigated such disks usi n g mid-IR high- 
resolution observations (Matsuura et al. 2006; Chesneau et al. 



120071: iLvkou et alj|201 UlLagadec et al.ll201 II) . 

Low-excitation rotational lines of carbon monoxide (CO) are 
particularly useful for probing the obscured cores at the center 
of evolved stars, including young PNe. They trace the mass dis- 
tribution, hindered by extinction in the optical, and provide di- 
rect information on the gas kinematics ([B uiarrabal et al. 2001, 
1998: IZweigle et al.lll997t iBachiller et aLl Tl988). Zweigle etal] 
(1 19971) mapped the 12 CO 7 = 2-1 line emission in M 2-9 with 
the Plateau de Bure Interferometer (PdBI). These maps revealed, 
with an angular resolution of 3 "x 5", a large expan ding equa- 
torial ring in the nebula waist. IZweigle et alJ (Q~997) estimated 
from 12 CO 7=1-0 line emission a mass of ~ 0.01 M . To ex- 
plore the equatorial mass distribution in M 2-9 in more detail, 
we carried out new PdBI high angular-resolution observations. 
We present in this paper the results of those observations, and 
analyze in detail the morphology and kinematics that lead to a 



new understanding of the nebula shaping and the central stellar 
system. 



2. Observations of 12 CO 7 =2-1 line and 1 mm 
continuum emission 

Observations of 12 CO 7 = 2-1 line emission (at 230.538 GHz) 
were performed with the Plateau de Bure interferometer in 
M 2-9. A track of observations was obtained in each Aq and Bq 
array configuration in February 2009. Channel maps were ob- 
tained with a synthetic beam (at half power) of 078 x 074 (PA. 
15°) in size. Image synthesis was performed with natural weight- 
ing and with the Hogbom cleaning algorithm. The results are 
presented in Fig.[T]with a channel spacing of 1 km s _1 , though 
a detailed data inspection was performed to spectral resolutions 
of 0.1 kms" 1 . In Figs. Q] and |2] the new M 2-9 CO data are 
superimposed ov er the optical ima ge obtained with the Hubble 
Space Telescope (IBalick et al.ll997l) . Astrometry was performed 
to properly overlay millimetric and optical data. 

Simultaneously, we observed a 1.8GHz bandwidth free of 
line emission, which was averaged to detect unresolved 1.3 mm 
continuum emission of 240 + 1 mJy in brightness. Note that this 
emission was subtracted in the interferometric visibilities be- 
fore imaging the CO line emission presented in this paper. The 
measured continuum flux is co mpatible with a previous analysi s 
of single-dish observations by ISanchez Contreras et al.l (Q998), 
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where the central compact continuum emission was deduced to 
originate from a region of warm dust close to the stellar system. 

In addition to the ring seen in previous studies (Zweigle et alJ 

1 1 997b . the new observations reveal a second tight inner ring in 
the waist of M 2-9. The new ring is almost three times smaller 
and its 12 CO 7 = 2-1 emission is about four times weaker than 
that of the outer ring. 

Fig- 121 presents in an inset the velocity-integrated 12 CO 
7 = 2-1 line emission synthesized with a circular beam of 074 
in size. Only the brightest CO line emission was included to bet- 
ter separate the closest parts of the two rings (northernmost and 
southernmost regions). As in Fig. [1] the inset is superimposed 
over the HST optical image. The continuum emission is repre- 
sented with a red dot. 

At first we deduce that the inclination angles of the two equa- 
torial rings coincide, and are perpendicular to the optical lobes 
(see further analysis in Sect. 13. lb . No emission is detected in 
the intermediate regions down to a level of 10 mJy beam -1 . The 
most extended molecular ring is detached by more than 1" from 
the inner one. Its outermost deprojected diameter is 7/2, and the 
innermost one 4'.'2. The easternmost part of the outer ring seems 
to reach regions outside the external optical lobes, perhaps with 
some flaring, and is thicker than the westernmost part. The inner 
ring is 2'.'4 in size, and its easternmost part also presents a mod- 
erate increase in brightness. From the spatial resolution and the 
signal-to-noise ratio, we consider our size estimates to be better 
than 072. 

Position-velocity diagrams of the 12 CO 7 = 2-1 line emis- 
sion in M 2-9 along the east-west and north-south directions are 
shown in Fig. [3] 



3. Data analysis 

Two coaxial rings or torus-like structures have been found in 
12 CO 7 = 2-1 line emission laying on the equatorial plane of the 
elongated M 2-9 nebula. 

3. 1 . Two outflows in M 2-9 

The relation of the two rings and the extended optical lobes is 
best shown in Fig. [2] This comparison suggests that the two 
rings are the equatorial counterparts of the two coaxial optical 
lobes (the outer one being more diffuse). The shaping of the rings 
(or torus-like structures) and the formation of the bipolar optical 
outflows seem hence to result from the same two events of mate- 
rial ejection. The low equatorial expansion velocities measured 
for the two rings (Sect. 13 .3b suggest that they likely originate 
from an increase in the mass loss when the optical lobes were 
formed. Note that by extrapolating the image of the innermost 
optical lobes, from the region where their image becomes diffuse 
(at ~ 2" from the equator) to the equatorial plane, we measured 
that the offset between the lobes and the rings in the equator must 
be smaller than 074 (the ring s being outward s; this relation was 
investigated in young PNe by Huggins 2003). 

Assuming that the rings are circular, we estimated the in- 
clination of their symmetry axes with respect to the plane of 
the sky from their projected and actual sizes (along the sym- 
metry and perpendicular axes). For the two rings we derived a 
similar inclination of ~ 19°, with an uncertainty of ± 2° due to 
ring asymmetries. The northern (southern) part of the ring is ap- 
proaching (receding from) us. Therefore, we see the south face 
of the ring. Similar inclinations have been derived fo r the sym- 
metry axis of the bipolar outflows seen in the optical ([Solf 2000; 
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East offset (arcseconds) 

Fig. 2. Velocity-integrated 12 CO 7 = 2-1 line emission is shown 
in the inset, superimposed over th e optical image o btained 
with the Hubble Space Telescope (Bali ck et all 1 19971) . Only 
the brightness distribution above 60 mJy beam -1 (which corre- 
sponds to ~ 8 a, where <x is the root-mean-square noise) was 
included to better separate the close projected image of the two 
rings. At the center the continuum emission is represented by 
a red dot, which traces the position of the stellar system. More 
details can be found in Sect. [2] 



ISchwarz et alj|199"7l) . This confirms that the rings are accurately 
placed in the nebula equator. 

3.2. The lifetime of the nebula: its formation by means of two 
ejections 

The position-velocity diagrams presented in Fig. [3] clearly show 
that both rings expand. By assuming that their expansion is 
isotropic, the systemic velocities corresponding to the large and 
small rings are 80.0 and 80.6 km s _1 . The systemic velocity as- 
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Fig. 3. Position-velocity diagrams of the 12 CO 7 = 2-1 line 
emission in M 2-9 along the east-west (on the top) and north- 
south (on the bottom) directions. The green (violet) dotted lines 
represent the kinematical and spatial center of the inner (outer) 
ring. Top: Integrated brightness from the whole nebula. Bottom: 
Brightness obtained within a slit along the projected symmetry 
axis. With a green (violet) solid line we plot the position-velocity 
gradient corresponding to the entire inner (outer) ring along the 
north-south direction. A white contour traces the emission in- 
tegrated over the complete nebula at low brightness level. Note 
that the (Vlsr) velocity axis shows the expanding gas speed pro- 
jected on the line of sight. A more detailed discussion can be 
found in Sect. 13.21 

sociated to the outer ring hence differs from that of the inner ring 
by ~ 0.6 km s _1 (+ 0.1 km s" 1 ). 

In Figs. [T]|3]we also see that the centers of both rings are 
shifted along the east-west axis. The inner-ring center seems to 
coincide with the position of the stellar system, as traced by the 
continuum emission, whereas the center of the larger ring is off- 
set eastwards by ~ 0f.'3. The uncertainty in this estimate mainly 
comes from departures from circular symmetry, since the east- 
ernmost part of the larger ring is thicker than the westernmost 
part. If we only consider the regions brighter than half the max- 
imum, we estimate that the measured offset is 0'.'34 + 0'.'07. 
Considering the weakest detected regions, the offset would be 
(X'3 1 + 0'.'09. Note that the uncertainty triggered by the interfer- 
ometer position precision is considerably smaller. 

These kinematical and spatial differences between the two 
rings strongly suggest that the velocity of the mass-losing star 
was different when it expelled the two rings. 

The ratios between the projected sizes and expansion veloci- 
ties presented in the position-velocity gradients along the north- 
south direction (in Fig. [3} provide the times at which the ejec- 
tions occurred (multiplied by the tangent of the inclination angle 



of the symmetry axis with respect to the sky plane). These dif- 
ferent gradients for the two rings (represented with solid lines 
in Fig. [3} would hence be another indication that the rings 
were formed at diffe rent epochs. For a distance of 650 pc (see 
ISchwarz et alJIl 99*71 and discussion in Sect. @), we deduce that 
the time elapsed between the ejections of the two molecular rings 
is 500 yr. This value cannot be much lower (by less than 10%) 
considering the uncertainties in the inclinations of the rings. The 
first ejection, which shaped the outer ring, would have happened 
about 1400 yr ago. The second ejection, responsible for the in- 
ner ring, would have occurred 900 yr ago. We estimate a maxi- 
mum uncertainty of 140 yr in the derived times considering all 
possible errors. As discussed there, sizes and times depend pro- 
portionally on the assumed distance; all other results presented 
in this paper are distance-independent. 

We have already argued that the optical and molecular data 
show different parts of the same nebula components, which were 
formed by the same (two) events of mass ejection. This implies 
that the two lobes seen in the optical are not coeval, but were 
formed by two ejection events separated by ~ 500 yr, contr ary to 
what was concluded in previous studies (Sm ith et alj 2005). 

The measured difference between the systemic velocities of 
the two rings is readily explained if the star responsible for the 
ejections changed its velocity. This can be justified if the star or- 
bits within a binary (or multiple) system. The very short time (~ 
500 yr) elapsed between the two ejections compared to the life of 
a solar-mass type star (~ 10 9 yr) strongly suggests that the same 
late AGB or post-AGB star was responsible for these events. 
We do not consider the remote possibility of having two distinct 
post-AGB stars to have formed the two outflows. Moreover, the 
time lag between the two outfl ows is ~ 5 times longer than the 
binary-orbit period of ~ 9 yr dDovle et al.ll2000t ICorradi et al] 
l20TTt lL 1V10 & Sokerl200ll) . This suggests that the ejections were 
triggered during a critical phase in the evolution of one star, such 
as ultimate and sudden increases in the mass-loss rate of the 
late AGB (or post-AGB) star. The interaction between the stars 
therefore would not have determined the time at which the axial 
ejections happened. This is supported by the kinematical offset 
detected between the two rings, which implies that the ejections 
occurred at two different relative positions within the binary or- 
bit (Sect. [33). 

Mass-loss variations have been observed in the envelopes 
around other AGB and post-AGB stars with differ ent time scales, 
from ~ 400 yr (of ~ 3 10" 6 M R yr 1 for R Cas; see lSchoierf2007t 
ICastro-Carrizo et alj 1201 Oh to ~ 2000 yr (of ~ 10" 4 M p yr"' 
for CRL618 and IRC +10 420; ISanchez Contreras et al l [2004: 
ICastro-Carrizo et al. I l2007h . some of them being also visible as 
arcs or rings in the circumstellar envelopes. In M2-9, from the 
estimate of ejection duration (Sect.|5]l, we derive a mass-loss rate 
of ~ 9 10~ 5 M yr -1 for the first molecular wind. The mass-loss 
rate of the second wind, ejected ~ 500 yr later, is likely lower. 
However, in spite of these evidences of mass-loss variations in 
the late AGB or early post-AGB phase, the current small sam- 
ple and significant differences prevent us from a more detailed 
analysis. 

3.3. The nebula kinematics 

The large and small rings expand with velocities of 7.8 and 3.9 
km s _1 (+ 0.1 km s _1 ). The expansion velocity of 7.8 km s _1 is 
moderate compared with that seen in most of the AGB outflows, 
but 3.9 km s is a very low velocity compared to AGB shell ve- 
locities. Recent studies of the inner eq uatorial regions of some 
post-AGB objects, such as CRL618 (San chez Contreras et al.l 
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• Star position in the binary orbit during the last ejection 

• Star position during the tirst ejection, assuming it was brief 
O Star position during the first ejection, assuming it was long 




Binary counterclockwise rotation 



Fig. 4. Sketch representing the position at which the ejections 
occurred in the binary orbit, according to the discussion in Sect. 
13.31 and the two scenarios presented in Sect. [5] regarding the 
first ejection. Note that at first we cannot discern the sense of 
the binary rotation, clock- or counterclock-wise (as seen from 
the Earth, from the south), and both possibilities are hence rep- 
resented in this figu re. A clockwise rotation is suggested by 
ICorradi et al.l d201 ll) and also by the position-velocity gradient 
detected in the southernmost part of the outer ring (see Sect. 

E2J. 



2004; lBuiarrabaletaIll2010h and M 1-92 dAlcolea et alJl2007h . 
show that equatorial outflows in the very late AGB (or early 
post-AGB) phase can expand at these very low velocities, which 
is also supported by our data. 

There is no sign of deceleration within the rings that could 
indicate some interaction with a slower previously expelled halo 
or surrounding interstellar medium: no deceleration position- 
velocity gradient is detected in the rings, the outermost layers 
of both rings present a good circular symmetry, and the mean 
expansion velocity is uniform in each ring. In addition, no CO 
emission was detected close to the rings or to the rest of nebula, 
which would indicate the presence of nearby interstellar molecu- 
lar gas. The good axis-symmetry of the optical nebula also sup- 
ports this interpretation. No circumstellar emission is detected 
either between the two rings and, moreover, the innermost ring 
is expanding much slower than the previous one, which renders 
any interaction between them impossible. We interpret the data 
to show that the rings expand freely, and therefore conclude that 
times derived from the measured velocities indeed correspond to 
times elapsed since the mass ejections. 

From the estimated binary period and assuming a stellar 
mass of the AGB star < 1 M Q , we deduce a binary distance < 
20 AU (corresponding to ~ 0703 at a distance of 650 pc), which 
is indiscernible in the present data. The measured position off- 
set between the ring centers of ~ 073 is therefore too large by a 
factor > 10 with respect to the orbit size. Therefore the relative 
motion of the rings, deduced from the difference in the velocity 
centroids (Sect. [3T2b . likely also triggers this offset. The millime- 
ter continuum emission should correspond to the position of the 
binary system (Sect. 0, which seems to coincide with the center 
of the inner ring. 

In Fig. [3] we see that the center of the inner ring remains 
very close to the binary position (from our perspective), a pos- 
sible displacement in the last 900 yr being smaller than our un- 
certainty in the determination of the ring center, which is esti- 
mated to be ~ 0.06" due to small ring inhomo geneities (slightl y 
larger than the position uncertainty, < 0.04". iReid et alJ [l988). 
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Fig. 5. East-west position-velocity diagrams (similar to those in 
Fig. [3} to compare the observations performed in 1997 (in red 
contours) with the new data observed in 2009 (in color scale 
and black contours). From left to right data are synthesized with 
different beams. Left: All contours and the image are synthe- 
sized with the beam corresponding to the observations obtained 
in 1997 (of 2738x1704 in size, PA 7°). Center: Contours are 
plotted with their respective original beams, of 1704 and 074 in 
size in the east- west direction for the 1997 and 2009 observa- 
tions, respectively. Right: All contours and the image are syn- 
thesized with the beam obtained in 2009. A detailed analysis is 
performed in Sect. [4] 

Accordingly we deduce that the inner ring should be moving in 
a direction closer to the line of sight than 12° (independently of 
the distance, by adopting a binary rotation speed of 1 km s _1 ). 
The last ejection happened therefore when the star was moving 
in a direction very close to the light of sight, in the easternmost 
or the westernmost part of the orbit. 

The relative offset between the two rings indicates that when 
the first ejection happened, the star was moving eastwards. Their 
systemic-velocity offset suggests that the last ejection occurred 
when the star was receding. However, at this point we cannot 
distinguish the sense of the binary rotation. If clockwise (as seen 
from the Earth, from the south), the first ejection would have 
happened in the semi-orbit that is farther from us and the last one 
in the westernmost part of the binary orbit. If counterclockwise, 
the first ejection would have occurred in the section of the orbit 
that is closer to us, and the last one in th e easternmost part. T hese 
two scenarios are represented in Fig. [4] Corradi et all (120111) de- 
duced that the observed east-to-west motion of the bright nebular 
knots takes place in the side of the bulbs facing the Earth. This 
would indicate that the binary presents a clockwise orbit (upper 
plot in Fig.|4|. Accordingly, we can determine the position of the 
stars when the ejections took place, as shown in Sect. 15.21 There, 
additional discussion supports the clockwise orbital scenario. 

4. Distance 

We adopted a distance of 650 pc, a value so far mostly con- 
sidered in the literature for M2-9 that reasonably agrees with 
the tentative ring ex pansion presented in this section. Recently 
ICorradi et al.l (1201 ll) proposed a farther distance of 1,300 pc. 

To investigate the ring-expansion proper motions, a compar- 
ison was performed between observations obtained in 1997 (un- 
published data by Neri et al.) and the new data presented in this 
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paper. In Fig.[5]we compare the east-west position-velocity dia- 
grams (equivalent to that in Fig. [3j for the two data sets, synthe- 
sized with different beams. A correlation analysis of the brightest 
parts in the outer ring clearly suggests expansion proper motion 
over the elapsed time. From the compared images in Fig. [5] we 
estimate a diameter increase that roughly ranges from 0.10" to 
0.24". We have performed a similar analysis by fitting the data 
in the wv-plane. The two data sets were fitted at the central chan- 
nel by two Gaussian functions of FWHP 1717 and 0781, for the 
eastern and western ring components, their positions being free 
fitting parameters. In this way we estimate that the ring diameter 
increased by 0725 ±0711. 

From the measured expansion velocity (7.8 kms -1 ) and 
the time elapsed between the two observations (11.5 yr), we 
calculated the distance to M 2-9 with the formula <i(pc) 
= 36.9/diameterJncrease("). Accordingly, from the proper- 
motion estimates we deduce that the distance to M2-9 ranges 
from 100 to 300 pc. By arbitrarily increasing our uncertainty to 
0720, the distance would range from 80 to 800 pc. Therefore, 
though we cannot perform a precise distance estimate, it seems 
improbable that the distance is larger than 1,000 pc. 

Because our detection of proper motion is only tentative, 
we have analyzed the uncertainty in the value derived by 
ICorradi et al.l (1201 1) . Distances shorter than 1,000 pc should not 
be discarded either from their analysis of the proper motions in 
reflected light, particularly considering the observational uncer- 
tainties and the inconsistencies in the measured offsets for the 
two lobes. 

The choice of distance has an impact on the main known 
nebula characteristics. The estimated sizes and times in particu- 
lar are proportio nal to the distance. For a distance twice larger 
(as propose d bv ICorradi et al. 201 lh . the total nebula extent (as 
mapped by ISchwarz et al.l 1997b would be ~ 2 10 18 cm, about 
three times larger than usual PNe sizes (e.g. iBuiarraba l et alj 
1988). M 2—9 would be one of the largest known PNe. Also, 
the two ejections would be separated by ~ 1,000 yr, the first 
one would have occurred 2,800 yr ago. We note that these 
times are extraordinary long from what we expect in the AGB- 
to-PN e evolution, which is estimated to last ~ 500 - 2,000 
yr dBuiarrabal et al.l fl988l) . This is particularly unexpected if 
we consider that M 2-9 is not thought to be a developed PN 
but an early PN dCastro-Carrizo et alJl200l[). based on the esti- 
mate of its effective stell ar temperature (ICalvet & Cohen|[T978l: 
LSwings& Andrillatlll979h . 

In view of the presented data comparison, the above discus- 
sion, and the uncertainties we have, we prefer to adopt a distance 
of 650 pc. This value is so far mostly adopted in the literature, 
and facilitates a comparison with previous results. 

However, we also note that all the main results obtained in 
this paper are independent of the adopted distance. Regarding 
the velocities, they are either directly measured or derived by the 
ratio of sizes and times. Therefore estimated speeds, and conse- 
quently masses, are distance-independent (Sect. |5). The incli- 
nation angle of the nebula symmetry axis is also derived from 
angular sizes (Sect. [XT)- For the estimate of the direction along 
which the inner ring moves (Sect. 13 .3b the distance dependence 
cancels out with the ejection lifetime. The angular size of the 
binary orbit would be smaller by adopting 1,300 pc, ~ 0.01", 
which would make it more difficult to resolve the binary with 
the current data (Sect. 13.31 ). again supporting our interpretation. 



5. Kinematics of the binary orbit 

From the properties of the two rings we can estimate the modu- 
lus of the orbital velocity of the star that ejected the rings, V. For 
this we assume that the orbit is circular, that V remained constant 
during the whole mass-loss phase, and that the same star ejected 
both rings. We note that no time variations have been reported in 
the rotation pattern seen in optical observations, and associated 
to the orbital rotation itself. Hence, no high eccentricities are 
expected in the bin ary orbit (as concluded by iDovle etall bOOO; 
ICorradi et al.l201 il) . which supports our main assumptions. Here 
we derive properties of the binary rotation from parameters di- 
rectly measured from the observations, and the mentioned as- 
sumptions. To start, note that the speed difference between the 
two rings allows us to give a lower limit for the orbital velocity 
of 0.31 km s (obtained by deprojecting the measured velocity 
difference, divided by two). 

Below, we refer to axis x as the axis perpendicular to the 
east-west direction within the equatorial plane, passing through 
the continuum position, its projection in the sky plane being par- 
allel to that of the nebular axis. Axis y is perpendicular to axis x 
within the equatorial plane, parallel to the east-west direction. 

5. 1 . Very short mass-loss event assumption 

In a first approach, we assume that the rings were ejected during 
phases much shorter than the orbital period. From the observa- 
tions, we deduce 1) A position offset between the ring centers of 
0734 (± 0707) in axis y, which results from the relative motion 
of the rings. 2) The second inner ring is moving in a direction 
closer than 12° to that of axis x, which is an upper limit. 3) The 
relative velocity between the two rings has a projection in the 
line of sight equal to 0.60 km s _1 (± 0.10 km s~'). From point 
(2) we can assume that the last ejection took place when the star 
was moving in the direction of axis x, so that the velocity of the 
inner ring was such that V' y ~ and V' x ~ V [Eq. El] (where V 
is the modulus of the binary rotation velocity). Then, by consid- 
ering item (1) the offset of 0734 in axis y is only due to the ve- 
locity of the outer ring, i.e. V° = 0734/r o = 0.74 km s -1 (± 0.23 
km s~') [Eq. E2], where T is the time from the first ejection (~ 
0.45"/km s _! from the gradient shown in Fig. [3] which multi- 
plied by the distance gives ~ 1400 yr). We hence obtain (V°) 2 = 
(V) 2 - (0.74 km s~') 2 [Eq. E3]. Finally, from point (3) we deduce 
that V° - V[ = 0.63 km s _1 (where 0.63 km s _1 is estimated by 
deprojection toward the equatorial plane of the measured 0.60 
km s" 1 ). By using Eq. El we obtain V° = V - 0.63 km s _1 [Eq. 
E4]. From the last two equations (E3 and E4) we deduce that V = 
0.75 km s _1 (± 0.50 km s~'). From V° and V we also deduce that 
the ejection of the outer ring happened from an orbital point at 
17° from the intersection with axis x. The velocity-vector direc- 
tion was hence close to that of axis y, though this angle estimate 
is uncertain within the southern semi-orbit. We finally note that 
assuming ejections shorter than they actually are would lead to 
underestimated velocities. 

We can also derive the systemic velocity of the binary 
stellar system in M 2-9 by correcting the measured sys- 
temic velocities of the rings by the derived stellar velocity at 
the time of the ejections. This is particularly straightforward 
from the inner-ring velocity. We obtain a stellar speed of 79.8 
km s" 1 (± 0.5 km s~'), which is slower and hence agrees bet- 
ter with the values recently s uggested from optical observations 
(iTorres-Peimbert et al. 2010). The fact that this value is lower 
than the outer-ring systemic velocity would indicate that the first 
ejection would have happened when the star was receding with 
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respect to us within the binary orbit, as represented in Fig. [4] 
though moving close to the east-west direction. 

5.2. Longer mass-loss event assumption 

Let us suppose that the first ejection happened over a longer 
time, during a significant fraction of the orbital period. From 
what is derived in Sect. 15.11 we can suppose that the oldest ring 
was ejected when the star was moving basically in the direction 
parallel to the east-west axis. The ejection happened when the 
star was describing an angle 2a in its orbit. The average velocity 
of the outer ejected ring in the y direction is then 

Vccscrdcr _ y s j„ 



2a 



[Eq. E5]. 



In previous estimates (Eq. E2) we derived that (V°) ~ 0.7 km s 1 
(+ 0.2 km s _1 ). Note that because of the longer mass-loss, we 
would expect an acceleration gradient in the ring along the per- 
pendicular axis (x) whose extreme values would be separated by 
AV X = 2Vsina [Eq. E6]. This would be because the gas expelled 
during the first part of the ejection phase would move faster in 
the expansion direction (x in our case) than the gas ejected subse- 
quently, as a consequence of the change in the ~V direction. This 
gradient should only be observed in the ring part corresponding 
to the same angle interval described by the star in its orbit during 
the ejection. Also, a compression effect should have happened at 
the opposite side of the ring (at 180°) due to the contrary relative 
movement between the different ring layers (see a descriptive 
sketch in Fig. [6]). 

We did indeed detect a small acceleration position-velocity 
gradient along the axis x in the southernmost part of the north- 
south diagram, shown in Fig. [3] the outermost layers of the outer 
ring expand faster than the innermost layers. Indeed, similar 
acceleration gradients are detected in the north-south position- 
velocity diagrams obtained in the central east-west 3"-size sec- 
tion. Considering the uncertainties we could have in our esti- 
mates of the gradient (regarding the beam size and the com- 
parison with the northern part) and also projection effects, we 
estimate that AV X = 2 km s _1 (+ 1 km s ). As predicted, we de- 
tected no sign of a velocity gradient in the northernmost side of 
the outer ring. From Eqs. E5 and E6 we derive a ~ 80°(± 60°), 
the ejection takes hence place during about half of an orbital pe- 
riod (+a) and V would be ~ 1km s + 1km s _1 , and therefore 
ranging from 0.3 (hard limit from the measured velocities; Sect. 
|5])to2kms" 1 . 



If our interpretation on the observed relative gradient in the 
southern part of the outermost ring is correct, this would be an 
additional indication that the ejection took place when the star 
was somewhere in the southernmost binary semi-orbit, suggest- 
ing a binary clockwise rotation (as shown in Fig. 0J. As de- 
scrib ed in Sect. | 3.3| thi s agrees with what was already suggested 
bv lCorradi et alJ(l201lh . 

We estimate this second scenario, based on the gradient de- 
tected in the ring southern part, to be more realistic, and hence to 
allow a better characterization of the first ejection and of the or- 
bital kinematics. We note, however, that a deeper analysis of this 
gradient is limited by the angular resolution in the north-south 
direction. 

From the deduced orbital velocity we estimate a systemic 
velocity for the binary stellar system of ~ 79.6 km s -1 , w hich 
also remains compatible with iTorres-Peimbert et alj d2010l see 
previous discussion in Sect. 15. II ). In addition, from the duration 
estimates we derive that the first ejection took place during ~ 40 
yr (+ 30 yr). For the inner ring, the determination of the ejection 
time is more difficult, but must anyhow be considerably shorter 
considering its narrowness. 

5.3. Implications for the companion mass 

We can estimate the mass of the companion from the derived 
orbital velocity. For this we assume again a circular binary or- 
bit and an orbital period, P, of ~ 90 yr dCorradi et alJ 1201 1: 
iLivio & Soke^l200ll iDovle et al-l feOOO). From the gravitational 
equations associated to the movement of two bodies we obtain 
that the modulus of the orbital velocity of the primary star (1: 
presumably the most massive and responsible for the mass ejec- 
tion) is in solar units: Vi(km s" 1 ) = 30 x (^f 3 [Eq. E7], 
where M m is the total mass and Mo that of the secondary. 

By introducing the derived V in Eq. E7, we find that the mass 
of the secondary must be very small. Indeed, for a mass of the 
primary (late AGB) star of ~ 1 M the mass of the companion is 
M2 ~ 0.1 - 0.2 M Q . If we adopt a different M tot by a factor 2, the 
mass of the companion would change by + 0. 1 M . Considering 
that the total estimated ci rcumste liar mass is t 0.1 M H (including 
ionized and atomic gas, ICastro-Carrizo et aUl2001l) . a mass for 
the primary much larger than 2 M does not seem very likely. 
However, note that if M to , would be ~ 4 M , the mass of the 
companion would be ~ 0.4 M Q . In all the cases the companion 
would be considerably smaller than the primary star. 



Ejected Ring 




Fig. 6. Sketch to visualize the origin of the velocity gradient ex- 
pected in one side of the ring, if its shaping took place for a 
fraction of the binary orbit period. According to the scenario de- 
scribed in Sect. 15.21 in the sourthern part of the ring the outer- 
most layers would expand faster than the innermost ones. 



6. Conclusions 

We presented subarcsecond-resolution mapping of 12 CO 
7 = 2-1 line emission in the proto-planetary nebula M 2-9. Our 
data show two independent expanding rings (or torus-like struc- 
tures) in the nebula equator. Continuum emission has been de- 
tected, likely tracing the position of the stellar system. 

The centers of the two rings are offset by 0.34", and their sys- 
temic velocities differ by 0.6 km s _1 . This is solid evidence that 
there is a binary system at the center of M 2-9. In this scenario, 
the two rings were very probably ejected by one of the stars, in a 
very late AGB or very early post-AGB phase. The different ve- 
locities of the star in its orbit when the material shaping the rings 
was expelled explain the measured position and velocity offsets 
between the rings. At the same time, these two mass-loss events 
likely gave rise to the two coaxial lobes seen in infrared and op- 
tical images The first ejection happened 1400 yr ago, the second 
500 yr later. We suggest that, after the material was expelled by 
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the star, the axial lobes were shaped, accelerated, by interaction 
with post-AGB fast and collimated jets. The interaction effects 
are still visible in the bright shock-excited knots detected in the 
visible but did not affect the dynamics of the equatorial rings, 
in agreement with current theoretical ideas of PN shaping (see 
Sect. 1). Although most of the material in the equatorial rings 
and the axial lobes was probably ejected by the primary (in a late 
AGB or early post-AGB phase), the structure and movements of 
the knots indicate that the colli mated exciting jets came from a 
compact, less bright secondary (Corra di et alJl201 lb . 

A distance of 650 pc was adopted in this paper, which agrees 
better with the tentative detection of proper motion we presented 
here. None of our main results (except for lifetimes) depend on 
the adopted distance. 

By analyzing the kinematics of the CO-emitting equatorial 
rings we estimated the positions at which the two ejections took 
place in the binary orbit. With some few simple assumptions, we 
derived an orbital velocity of ~ 1 km s _1 . From the equations of 
the orbital movement of two objects, we obtained a very small 
mass for the secondary, J 0.2 M Q (+ 0. 1 M Q , by assuming that 
the mass of the primary star is between 0.5 and 2 M ). We sug- 
gested that the binary system is composed of a red giant and a 
dwarf. We estimated that they would be separated by 20 ± 5 AU 
(for the interval of total masses considered above), which is too 
large for the stars to have undergone a com mon envelope phase , 
or for being interacting symbiotic systems fMikolaiewska 2007; 
iPodsiadlowski & Mohamedll2007l) . 

The strongly axisymmetric nebula and very collimated jets 
around the binary system convincingly suggest that the multiple 
nature of the stellar component played a relevant role in the ejec- 
tion and shaping of M2-9. The significance of this low-mass 
companion to shape the extended and elongated bipolar nebula 
seen in M 2-9 has important consequences on our understand- 
ing of the mechanisms that drive the shaping of asymmetrical 
planetary nebulae, since the existence of companions is increas- 
ingly more probable when the requirements for their masses are 
relaxed. 
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